library(foreign)
library(lmtest)
library(sandwich)
new.data <-read.dta("C:\\Documents and Settings\\anhill\\My Documents\\Gov 2001\\Replication\\New Data.dta")
ls()
head(new.data)
names(new.data)
attach(new.data)

new.lm21 <-lm(z_femlabor06 ~ z_loggdpcap06 + z_loggdpcap206 + z_workingage06
+ z_me_nafr + z_commie)
new.lm21
summary(new.lm21)
coeftest(new.lm21, vcov=vcovHC(new.lm21, type="HC1"))

new.lm22 <-lm(z_femlabor06 ~ z_loggdpcap06 + z_loggdpcap206 + z_workingage06
+ z_me_nafr + z_commie + z_islam)
new.lm22
summary(new.lm22)
coeftest(new.lm22, vcov=vcovHC(new.lm22, type="HC1"))

new.lm23 <-lm(z_femlabor06 ~ z_loggdpcap06 + z_loggdpcap206 + z_workingage06
+ z_me_nafr + z_commie + z_oilrentspercap06)
new.lm23
summary(new.lm23)
coeftest(new.lm23, vcov=vcovHC(new.lm23, type="HC1"))

new.lm24 <-lm(z_femlabor06 ~ z_loggdpcap06 + z_loggdpcap206 + z_workingage06
+ z_me_nafr + z_commie + z_islam + z_oilrentspercap06)
new.lm24
summary(new.lm24)
coeftest(new.lm24, vcov=vcovHC(new.lm24, type="HC1"))


### TABLE 4 ###

new.lm41 <-lm(z_femaleseats06 ~ z_loggdpcap06 + z_me_nafr)
new.lm41
summary(new.lm41)
coeftest(new.lm41, vcov=vcovHC(new.lm41, type="HC1"))

new.lm42 <-lm(z_femaleseats06 ~ z_loggdpcap06 + z_me_nafr + z_islam)
new.lm42
summary(new.lm42)
coeftest(new.lm42, vcov=vcovHC(new.lm42, type="HC1"))

new.lm43 <-lm(z_femaleseats06 ~ z_loggdpcap06 + z_me_nafr + z_oilrentspercap06)
new.lm43
summary(new.lm43)
coeftest(new.lm43, vcov=vcovHC(new.lm43, type="HC1"))

new.lm44 <-lm(z_femaleseats06 ~ z_loggdpcap06 + z_me_nafr + 
z_islam + z_oilrentspercap06)
new.lm44
summary(new.lm44)
coeftest(new.lm44, vcov=vcovHC(new.lm44, type="HC1"))

new.lm45 <-lm(z_femaleseats06 ~ z_loggdpcap06 + z_me_nafr + 
z_islam + z_oilrentspercap06 + z_femlabor06)
new.lm45
summary(new.lm45)
coeftest(new.lm45, vcov=vcovHC(new.lm45, type="HC1"))

new.lm46 <-lm(z_femaleseats06 ~ z_loggdpcap06 + z_me_nafr + 
z_islam + z_oilrentspercap06 + z_polity)
new.lm46
summary(new.lm46)
coeftest(new.lm46, vcov=vcovHC(new.lm46, type="HC1"))

new.lm47 <-lm(z_femaleseats06 ~ z_loggdpcap06 + z_me_nafr + 
z_islam + z_oilrentspercap06 + z_polity + z_pr_mod)
new.lm47
summary(new.lm47)
coeftest(new.lm47, vcov=vcovHC(new.lm47, type="HC1"))

new.lm48 <-lm(z_femaleseats06 ~ z_loggdpcap06 + z_me_nafr + 
z_islam + z_oilrentspercap06 + z_polity + z_pr_mod 
+ z_closed_list)
new.lm48
summary(new.lm48)
coeftest(new.lm48, vcov=vcovHC(new.lm48, type="HC1"))


### Table 5 ###

new.lm51 <-lm(z_femalemin05 ~ z_loggdpcap06 + z_me_nafr)
new.lm51
summary(new.lm51)
coeftest(new.lm51, vcov=vcovHC(new.lm51, type="HC1"))

new.lm52 <-lm(z_femalemin05 ~ z_loggdpcap06 + z_me_nafr + z_islam)
new.lm52
summary(new.lm52)
coeftest(new.lm52, vcov=vcovHC(new.lm52, type="HC1"))

new.lm53 <-lm(z_femalemin05 ~ z_loggdpcap06 + z_me_nafr + z_oilrentspercap06)
new.lm53
summary(new.lm53)
coeftest(new.lm53, vcov=vcovHC(new.lm53, type="HC1"))

new.lm54 <-lm(z_femalemin05 ~ z_loggdpcap06 + z_me_nafr + z_islam
+ z_oilrentspercap06)
new.lm54
summary(new.lm54) 
coeftest(new.lm54, vcov=vcovHC(new.lm54, type="HC1"))

new.lm55 <-lm(z_femalemin05 ~ z_loggdpcap06 + z_me_nafr + z_islam
+ z_oilrentspercap06 + z_femlabor06)
new.lm55
summary(new.lm55) 
coeftest(new.lm55, vcov=vcovHC(new.lm55, type="HC1"))